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Abstract 

We discuss CLT for the global and local linear statistics of random 
matrices from classical compact groups. The main part of our proofs 
are certain combinatorial identities much in the spirit of works by Kac 
and Spohn. 

1 Introduction 

Let M be a unitary matrix chosen at random with respect to the Haar 
measure on the unitary group U{n). We denote the eigenvalues of M by 
{exp(i ■ where — vr < 61,62, ... ,0n < ti- The joint distribution of the 

eigenvalues (called the Weyl measure) is absolutely continuous with respect 
to the Lebegue measure 11^=1 "^^i n-dimensional tori and its density 
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is given by 

Puin){eu...,9r.)^ } ■ n |exp(z-^,)-exp(z-^fe)|2 (1.1) 

^ ' l<j<k<n 

(see [We]). Throughout the paper we will be interested in the global and 
local linear statistics 

n 

'5n(/) = $^/(^,), (1.2) 



Sn{g{Ln-)) = Y.g{Ln-e^), (1.3) 
r n 

L„ ^ OO, )• 0. 

n 

The optimal conditions on /, g for our purposes are 

oo 

E l/»MA:|<oo, (1.4) 

fe=— oo 



\g{t)\^ ■\t\dt < oo, (1.5) 



where 



fix) = />) • ^''^ 

fe=— oo 

1 r°° 

g{x) = ^ • / ^(t) • e^*^dt 

V^TT J-oo 



However in order to simphfy the exposition we will always assume that / has 
a continuous derivative on a unit circle ( / G C^(S'^) ) and is a Schwartz 
function ( g E /(J(M^)) ). 

Let us denote by the mathematical expectation with respect to Haar 
measure. Wc start with the formulation of the result which is essentially due 
to C. Andreief ([A], for a modern day reference see [TWjand also [Dy]. 
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Proposition 

E^expitSnif)) - det(/rf+ (e*^ - 1)K,) = det(/d+ (e*^ - (1.6) 

where (e*-^ — 1) is a multiplicaiton operator and Kn,Qn '■ L^{S^) — > L'^{S^) 
are the integral operators with the kernels 

1 sin (^(x — y)) , , 

Kn{x,y)^- )')''\ (1.7) 



n-1 ^ ^ 

g„(x, y) = 5] ^e^'^'^e-'^y (1.8) 

Rerricirk 1. Kn,Qn are unitary equivalent to each other and are the oper- 
ators of a finite rank. In particular, Qn is just a projection operator on the 
first n harmonic functions of the unit circle. 

One of the ingredients of the proof of the proposition is the following 
chain of the equalities 

Pu{n){Ou . . . ,^n) = ^ ■ det(e^-(^-^>^'=)i<,-fc<„ • det(e-^-(^-^>^'=)i<,- fc<„ 

= -^det(g„(6',-,^fe) ) ^ gN 



1 , 

— det 
nl \ 



l<jr',fe<ra 



Remark 1 allows us to rewrite the Frcdholm determinants in (1.6) as the 
Toeplitz determinant with the symbol exp(t • /(•)): 



Enexp(^tf2f{e,)^ =D„_i(^exp(t-/)^ 



det( — / exp(t/(x)) • exp(i(j — k)x)dx 1 

P-^ Jo y 'i-<j,k<n 

(1.10) 
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The asymptotics of (1.10) for large n is given by the Strong Szego Limit 
Theorem: 

D„_i^exp(t • /) 

. +00 . (1.11) 

expUnm + -t'Y,\k\\f{k)\' + 0{l)) 
^ —00 ' 

(see [Sz] and [K], [H], [De], [F-H], [G-I], [Widl], [Wid2], [McC-W], [Ba-W], 
[Jol], [Bo], [Bo-S], [Me], [So2], [Wie], [D] for further developments.) 

In probabilistic terms (1.11) claims that ESn{f) = ■ j\ f{0)de + 0(1) 
(actually the remainder term is zero), and the centralized random variable 
YTj=i fi^j) = Y^]=i fi^j) converges in distribution to the normal law 

Our first goal is to establish a similar result for the local linear statistics. 

Theorem 1. Let g G J(M^), L„ +00, ^ ^ 0- Then EnY.%19i.Ln ■ 
^i) ~ 2t^l ' I'^oodi^)^^' ^''^'^ centralized random variable Y^iidi^n ■ 
9j) — EYTj=ig{.LnOj) converges in distribution to the normal law N{0, ^ ■ 

Clgmtldt). 

We give a combinatorial proof which holds both in the local and global 
cases. In some sense our approach is close to the heuristic arguments in [I-D] . 
We start with 
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Lemma 1. Let Q,n(/) be the i-th cumulant of Sn{f)- Then 



\m— 1 



m 

ki+...+k(=0 m=l 



£1 / ti h+(^2 

I /7 _ /7 ^' ^ i=l i=l 



ll + ... + im=^, 

i=l 



ei>i,... 


/m>l 




-1 ^1 


E 


— max(0, X^(- 


1=1 


1=1 


il+...+tm. 




E 

i-l 


(-^i))^ 1 < const^ 



\h\\f{h)\-...-\m\ 

(1.12) 



ki + ... + ki = Q 
|fei|+...+|fe^|>n 



Remcirk 2 One can see that for sufficiently smootfi / tfie r.fi.s. of (1-12) 
goes to zero as n — > oo. 

Remark 3 An analogous result to lemma 1 was established in [Spo] for the 
determinantal random point field with the sine kernel (see also Remark 4 
below). 

The proof of Lemma 1 will be given in §2. At this state we observe that 
it implies 

Lemma 2 The limit of C£^n{f),i > 1 exists as n — > oo and is equal to 
^ki+...+ke=o f(ki) ■ ■ ■ ■ ■ hh) ■ (G(ki, . . . ,ke)+G(-ki,... ,-ke)), where G is 
the piece-wise linear continuous function defined by 

1)"^ 1 



ae5, m=l + . . . + £^ = 

^1>1,.../^>1 (1.13) 

max I 0, X kain)., ^ K{i), . . . , ^ K{i) 

i=l 1=1 1=1 
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Proof of Lemma 2 After opening the brackets in (1.12) we observe that 
the coefficient in front of n is equal to 

y y (-i)"- « = Ji- « = i (1.14) 

£i>l,i=l,... ,m 

Indeed, the generating function of these coefficients is equal to 

logfl + (e^-l)^ 

□ 

Now CLT for YTj=i fi^j) follows from 

Main Combinatorial Lemma 

Let ki, . . . ,kj^ be arbitrary real numbers such that their sum equals zero. 
Let G{ki, . . . ,ki) be defined as in (1.13). Then 



G{ki, ... ,ki) 



= Ifel ifi = 2 

ife>2' 

We will prove the lemma in §3. 

Remark 4 A similar combinatorial lemma was stated by Spohn in [Spo] . He 
studied a time-dependent motion of a system of infinite number of particles 
governed by the equations 

i^j ~ 

where {6j(t)}^^^-independent standard brownian motions, and the initial 
distribution of particles is given by determinantal random point field with 
the sine kernel However, no correct proof of the combinatorial 

TT{x—y) ' '^^ 

result was given there. For completeness wc give a proof of Spohn's lemma 
independently from the proof of our Main Combinatorial Lemma in §3. 

Assuming the combinatorial part is done we can quickly finish the proof of 
Theorem 1. The formula for the mathematical expectation is trivial. Rewrit- 
ing (1.12) for the higher cumulants of Yl]=i9i^n • ^j) we see that the limit 
of the £-th cumulant is given by 

(27r)-i ■ j g{ti) ■ . . . ■ g{te) ■ (^G(ti, . . . ,te) + G{-ti, ... , -te)^ dt^... die 



6 



where the integral is over the hyperplane ti + . . . + 1£ — 0. 
Theorem 1 is proven. 



□ 



Remark 5 Our method also gives an elementary combinatorial proof of 
Szego theorem ((1.11)) for / G C^{S^) and sufficiently small complex t. It is 
different from the one suggested by Kac in [K] where the Taylor expansion of 
Dn{l — tg) as a function of t was calculated and then a so- called Kac-Spitzer 
combinatorial lemma was employed to confirm (1.11). 

Remcirk 6 Results similar to Theorem 1 have been established for other 
random matrix models in [Spo], [Jo3], [KKP], [Ba], [B-F], [SSol], [SSo2], 
[BM-K]. 

The rest of the paper is organized as follows. We prove Lemma 1 in §2 
and Main Combinatorial Lemma in §3. The result analogous to Theorem 1 
for orthogonal and symplcctic groups is established in §4. 

The author would hke to thank Ya. Sinai, P.Diaconis, K. Johanson and 
A. Khorunzhy for useful discussions. The work was partially supported by 
the Euler stipend from the German Mathematical Society. 

2 Proof of Lemma 1 

We start with calculating the moments of Sn{f). Le us remember that k- 
point correlation function of the eigenvalues of random unitary matrix is 
given by 

n\ f 



where the indices ii, . . . ^in range independently from 1 to n, and in par- 
ticular can coincide. Let M. = {Mi,... , M^} be a partition of the set 
{1,2,... , N} into subsets determined by coinciding indices among ii, . . . , ijv : 




(2.1) 



The A^-th moment of Sn{f) is equal to 
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Ml = {jl'\...,£^},...,Mr = {A'\...,jt^}M=iMi = {1,2,. ..,7V}, 
Si = |Mi|,z = 1, . . .r. Then 

^n(^n(/)) = E E (2-2) 

^ ^ overall ^it^^zt^.-t^^ 

partitions M 

Let US consider a typical term in (2.2) corresponding to a partition 



En J2 r(^^J-----r(^J = 

/ rH^l) ■ ■■■■ r^'M ■ PnA^l^ ■ ■ ■ ,Xr)dXi ...dXr 

By definition of the determinant and (2.1) 

r 

pn,r{^l,--- ,Xr)^ ^ (-1)'" (5„(Xi, X^(j)). 

Writing the permutation cr e S'r as a product of cychc permutations we have 

over partitions 

Pa X V y 



Fa X 
nil nvoMo .7—1 



over all cyclic 
permutations of /fa 

where {1, . . . ,r} = UfKa.Ka = {t^"\ . . . ,4°^}, a = 1, . . . ,q,Pa = \Ka\. 
Substituting (2.4) into (2.3) we arrive at the expression that has the following 
form : 

over partitions over partitions 

M={Mi,...,Mr} of {l,...,Af} lC={Ki,...,Kq} of {1,... ,r} 

To interchange the order of summation we construct a new partition V = 
{Pi, ... , Pg} of {1, 2,...N} as follows: Pi = Ujg^^.Mj, i = 1,... ,q. Then 
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{Mj}jQKi gives a partition of Pj that we denote by Vi- We have 

En{sM)r= E (n( E 

over partitions ^ over partitions 

V={Pu-,Pq} o{{l,...,N} Vi of Pi:n={P,,i,...,P,,t 



over cyclic 
permutations <T&St^ 

Yl Qn{Xj, Xa(j))dXi . . . dxt^) j . 

j=l / 

We remind that the moments are expressed in terms of cumulants 

mjv= ^ C\p,\ ■ . . . ■ C\Pi^\. 

over partitions 
•p={Pi,...,Pfc} 

Comparing the last formula with (2.5) we arrive at 

partitions 
r={Ri,...,Rm} of {!,... 



cyclic permutations-?—^ 

aeSn 

^ over ordered collections 
{ei,...,em).ET^i='^A>^ 

m 

Qn{Xj, Xj^i)dXi . . . dXjn 



= E E 



m=l 



^j>l,i=l,... ,m 
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m- JT' 



JJ^ Qni^jj Xj^i)dXi . . . dXm- 

Since Qn{x,y) = ^"Zq e^*-'^^"^) we can rewrite (2.7) as 



f^ix,)-...-f-ix^)- 

(2.7) 



CiAf) = E E 



^ / ^m— 1 



€l+...+^m=A 4>1 

n— 1 n— 1 



E • • • E /''(-^- + • + ^2) • . . . • /^-(-Srn-1 + Sj. 



S1=0 Sm=0 



Writing down the Fourier coefficients of the powers of / as the convolutions 
of the Fourier coefficients of / 

(fci, . . . ,kij : 

kl+...+k(-^=Sl—Sm 

F^(-si + S2)^ E • • • • • • • • 

{hi+i, ■ ■ ■ ,^2) ■ 

f"^{-Sm-i + Sm)^ E f{ke^-i+i) ■ ■■■■ hkej, 

we obtain 

e 

/ 1 \m—l 

CeAf)- E />i)-----/>.)-EH^ ^ 

fei+...+fe,=0 m=l (£1,... : 

#{M:0<M<n-l,0<ii+ ^^-^^ 



^1 • • • • ^m- 

£1 ei + .-.+im-l 

'^ki<n-l,... ,0<u+ E 
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The last factor in (2.8) is equal to 

(£i ei+...+em-i \ / ^1 ei+...+em-i 

0,^ki,..., ^ kj j - max I 0, y^^{-ki), ... , ^ {-ki 

(2.9) 

if the expression in (2.9) is nonnegative or zero otherwise. 

Lemma 1 is proven. □ 



3 Proof of the Main Combinatorial Lemma 

First we show that G{ki, ... , A;^) is a linear combination of terms \kij^ + . . .+ 
ki^ I . Then we compute the coefficient in front of every such term and show 
it to be equal to zero. 

Assume £ > 2. Consider a partition V = {Pi,... , Pm} of the set 
{1,2,... Let us denote Vi = X^jePi • • • ' ~ SjeP^ ^j- The ex- 

pression for G can be transformed into 



m=l V={Pu...,Pm} r&Sm ^ ) 

max(0, Vr{l),Vr(X) + Vr{2), • • • , ^'r(l) + Vr{2) + • ■ ■ + t'r(m-l)- 

In [R-S] Rudnick and Sarnak, following the ideas of [K] and [Spi] (see also 
[B], [An] ) , used the following identity for the set of real numbers Vi, . . .Vm 
with zero sum: 

- ^ max(0, Vril), Vr{l) + Vr{2), ■■■ , Vr{l) + Vr(2) + ■■■ ^^r(m-l) 



m 

reS, 



(3.2) 

> I I r I — \ \i\ 'rri, — \ r \ — i 1 1 • i > 'iii>\ 

Fc{l,...m}, «eF 



\ (|F|-l)!(m-|F|-l)!•|;^^;,| 
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The last formula gives us 



m=l V={Pi,...,Pm} Fc{l,...m}, 

F,FC=illi 



(3.3) 



(1^1 -1)! 



(m- |F| - 1)! 



Le us denote by A the subset Uj^pPj of {1,2,... , i}. Then {Pj}j(zp defines 
a partition of A, and {-PjjjeFC a partition of A*" = {1,2,... A. 

Wc change now the order of summation in (3.3): first we sum over all 
nonempty subsets A of {1,2, .. . ,i} and then over all partitions of A and 
A^: 



Ac {!,... ,£}, over partitions 

A,AC^ii U={Ui,... ,Ur} of A 



E 



{\U'\-1)\ 



over partitions 
W of AC 



ieA 



(3.4) 



Finally we note that 



1^1 



E (-i)i"i-^-(iwi-i)! = E E (-1) 

ELi«i=l^l>*i>i 

|A|! 1 



W={f7i,...f7^} 

(r-1)! 



the expression we already considered in (1.14). Indeed, there are exactly 
^J^^'',^ I ■ ^ different partitions of A such that {ti, . . . = {|f/i|, . . . , \Ur\}. 
If jyl| > 2 this sum is zero. If \A\ = 1, then {A'^l = i - \A\ > 2 and the 
second factor in (3.4) equals zero by the same argument. □ 
Now we turn to a combinatorial lemma first formulated in [Spo]. Let 
us denote by a = {ai, . . . , a^), /3 = (/3i, . . . , P^) vectors with entries aj £ 
{0, 1}. We consider a lexicographic order on the set of such vectors: a < {3 
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iS aj < l3j, j = !,...,£ and at least for one jo o^jo < Ao- Following 
[Spo] we call such nonzero vectors branches and a set T of ordered branches 
T = a(^) < a(2) < ... < a^^\ \T\ = m < £, a. tree. 

We denote by T{i) the set of all trees formed by a £- dimensional vectors 
(branches). A combinatorial sum in question is 

U{ki, ... ,ki)^ (-1)1^1-^ • max(0, a ■ k\a e T) (3.5) 

Here we used the notation a ■ k — o-j ■ kj. We call max(0, a ■ k\a E T) 
the maximum of the tree T. For a warm-up we prove 

Proposition 1 

U{ki,... ,ke) + U{-ki,... ,-k(-) = G{ki, . . . ,ki,ki+i) 
+ G{-ki, ... , -ki, -ke+i), 

where kj^+i — —ki — k2 — . . . — k^. 

Remark 7 Once the proposition is proven we see of course that U{ki, . . . ,k(i)+ 

U{—ki, . . . , —kf) is zero for i > 2. 

Proof of Proposition 1 In the above notations 

(-1)''^'"^ 

G{ki, . . . ,ke,ke+i) = — max(0, a ■ /c'|q; e T), 

TgT(^+1) ' ' 

where k' = (ki, . . . , k^, A;^+i), and the sum is over all trees T G T(i + 
1) such that the largest branch of T, a^^^^^ is less than D = (1,1,... ,1). 
Similarly, we can write U{ki, . . . ,ke) = X^xer(£+i)(~-'-)'"^'^^ ■max(0, a-k'\a e 
T), where the sum is over the trees T e T{£ + 1) such that the {£ + 
coordinate of a*^'^') is zero. Wc define a "rotation" on the set of all 
trees such that a^^D ^ D : W^((«W, a^^), . . . ,a^\^^^)) = («(') - _ 
a^^\ . . . , a^l-^') — a^^\ D — a^^^). Since Yl^j=i kj = 0, we observe that 

max(0,Q; • k'\a e T) +max(0,Q; • {-k')\a e T) = 

max( 0,a- A;'|a G W{T) \ + max(0, a • {-k')\a G W{T) 




6) 
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The last equality implies 

U{ki, ... ,ke) + U{-ki, ... , -ke) = {-1^^^-^ ■ (max(0, a ■ k'\a e T) 

TeT{e+i) 

+ max{0, a ■ {-k')\a e T)) 

1 1^1-^ 

(max(0,a- G WP{T)) 

+ max(0, a ■ {~k')\a G WP{T)) 

(-l)\T\-l 

= (max(0,a • A;'|q; e T) 

rer(£+i) 

+ a-{-k')\a e T)) 

= . . . , ke+i) + G{-ki, ... , -ke+i) 

Here we used that for any T' with a^'-^'l) ^ D there exist a unique T with 
ajljj) = and < p < |T| such that T' = W^^'(r). □ 

Proposition 2 

;7(A;i,... ,A;f) = 0if£>2. (3.7) 

□ 

We proceed by induction. 

It is easy to check the case i — 2. Let us assume that the proposition 
is true for some i > 2. Consider U{ko, ki, . . . , ke). Since f/ is a symmetric 
function we may assume ko < ki < . . . k^. The continuity of U implies that it 
is enough to check (3.7) for nondegenerate vectors (/cq, ki, . . . ,ki). Therefore 
we may assume that the coordinates ki, . . .ki are linearly indepdendent over 
the integers. Fix such ki, . . .ki and consider [/ as a piecewise linear function 
oi y = ko, U{y, k) = U{y, ki, /c2, • • • , kg). Our first claim is that U{y, k) is 
zero for all negative y. To show this we write 

Uiy,k)^ J2 {-lf^-'m^iO,a-iy,k)\aeT) 

TeT{£+l) 

= E + E + E 

TeT{e+l): TeT(£+l): TGT{e+l): 
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We denote the three subsums by ?7i, C/2, The first subsum is equal to 



and by the induction assumptions both are zero. Now we spht the third 
subsum in two. Consider the smallest branch a G T such that the first 
coordinate of a is 1, denote this branch by a' and denote the preceding 
(may be empty) branch by a" . We write — U^^i + where in C/31 
the summation is over T e T{i+ 1), such that cci^'^ = 1, a^'^'"^^ = 1 and 
a' — a" > (1,0,... ,0), and in 1/3^2 the summation is over all other trees from 
U3. We establish a one-to-one correspondence between C/31 and C/3^2: for any 
tree Ti with a' — a" > (1,0,... ,0) we construct T2 = {a*-^-*, . . . , a", a" + 
(1,0,... , 0), a', . • • ) a'^'}- Clearly, IT2I = |Ti| + 1, therefore 



.niax(0,a- (2/,A;)|a G Ti) = -(-1)1^^1-^ • max(0, a • (y, A;)|a e T2), 



and C/31 and C/3 2 cancel each other. 

Now we assume that y is nonnegative and < y < ki < k2 < ■ ■ ■ < ki. As 
we already noted U{y,ki, . . . , kg) is a piecewise linear continuous function. 
We claim that it can change its slope only a,t y — 0. Indeed, U{y, fci, . . . ,ke) 
can change its slope only at the points of degeneracy of {y, ki, . . . , kg), where 
Oq ■ y + a. ■ k = ■ y + ■ k and the coordinates of (ag, a), (a^, a') take 
values zero and one. Because is a non-degenerate vector we must have 
y + a-k — a'-k{oT:a-k = y + a'k) . Since the tree T contains both branches 
(1,0;) and (0,0;') only if a' < a, the only solution for nonnegative vector 
(y, k) must he y = 0, a' = a. We will finish the proof of the proposition if 
wc show that U{y, k) = for sufficiently small positive y. We again write 
U = Ui + U2 + Us as before. Then L^i = by inductive assumption and Us 
is zero for sufficiently small y (C/3^1 and C/3^2 still cancel each other). We can 
write the second subsum C/2 as 



max(0, a ■ k\a e T), 



TeT{e) 



the second 




TeT{tj 




(3.8) 
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( the last sum includes empty tree). The first term in (3.8) is zero by inductive 
assumption and the second is also zero since 

V (-1)1^1= V (-1)-^ (^ + 1)! _o 

ti>l 

Proposition 2 is proven. □ 



4 Orthogonal and symplectic groups. 

We start with the orthogonal case. The eigenvalues of matrix M e S0(2n) 
can be arranged in pairs 

exp(i6'i),exp(-i6'i), . . . , exp(i6'n), exp(-z6'„), < 6*1, 6*2, . . . , 6'„ < tt. 

Consider the normalized Haar measure on S0(2n). The probability distri- 
bution of the eigenvalues is defined by its density (see [We]): 

Pso2n{ei,... ,e^) = 2-(^- J] (2cos^,-2cos^,)2 (4.1) 
The /c-point correlation functions are given by (see [Sol] ) 

Pn,k{Ol: • • • ,^fe) = detfx^n-ll^.^^)) (4-2) 

V / l<i.i<n 



where 



^2^-1(2;, y) = K2n-l{,X, y) + K2n-l{,X, ~y) = 

1 /sin( (^"-^f-^) ) sin( (^-y+^) )^^ (4.3) 
^' [ sin(^) + sin(^) 

In [D-S] and [Jo2] Diaconis-Shahshahani and Johansson studied asymptotic 
properties of linear statistics Yl^=i fi^j) where for simplicity we may as- 
sume that / is real even trigonometric polynomial, f{9) = Yl^=i '^ki^^''^ + 
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^ ^^^), ttk = /(k), k = 1,2, . . . ,m. As before we denote the linear statistics 
by Snif). Then = Trace {Y^'^=i ^kM''). It was shown that 



E2nex.plt-J2fi^j) 



' (4.4) 
exp ( ^2 ^(1 + (-!)')/» + ^ E ^/»' + 0(1) 

k=\ fe=l 

which imphes the convergence in distribution of ^YTi=\ fi^j) normal 
law 



(-. m m 
fe=i fe=i 



(Actually (4.4) holds under much weaker conditions — it is enough to assume 
/eCi+°([0,7r]),a>0). 

Remark 8 Similarly to the unitary case (4.4) is equivalent to the large 
n asymptotics result for some determinants, this time Hankel determinants 
(see [Jo2], [Jol]). 

Our combinatorial approach allows to prove CLT for all / G C^([0, tt]) as 
well as to study the local linear statistics g{Ln ■ {Oj — 0)),0 < 9 < Tr.In 
particular we establish 

Theorem 2 Let g be a Schwartz function, Ln — > +oo, — ^ and < ^ < tt. 
Then 'YTj=i di^n ■ {Oj — = ■ g{x)dx + 0(1), and the centralized 
random variable Xlj=i Qi^n ■ (% — 9)) — X]j=i 9'(-^n ■ (^j — ^)) converges 
in distribution to the normal law N{0, ^ |^(t)p|t|(it). 
Theorem 2 also holds for S0(2n + 1) and Sp(n). 

Let M G S0(2n + 1). Then one of the eigenvalues of M is 1 and the 
other 2n eigenvalues can be arranged in pairs as before. The density of the 
eigenvalues is equal to 

^SO(2n+l)(^l,... ,^n)= (-) " J] (2 COS - 2 COS ^,)2 • f[ Siu^ 

l<i<j<n i=l ^ ^ 

(4.5) 

The formula for the A;-point correlation function is 

PnA^i, ■■■,0k) = det {K^^{9„ 9,))^ , (4.6) 
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where 

K2nix, y) = K2nix, y) - Kinix, -y) = 

1 I sin {n{x — y)) sm{n{x + y))\ (4.7) 
2^[~^^^) sin (2^) 

The analogue of (4.4) reads 



( m 2 ™ _ ^ 

^2 ^ + ^-^)') ^^^^ + 2 5^ + 
fc=i fe=i / 



(4.8) 



In the symplectic case M e Sp(n) the 2n eigenvalues again can be arranged 
in pairs 

exp(z • 6'i), exp(-i • 6'i), . . . , exp(i • 6'„), exp(-i ■ 9n),0 < 61,62, ■ ■ ■ ,6n < 7^, 
their density is equal to 

Psp(n)(6^,... ,6n)= i-j ■ n (2cos^,-2cos^,)2-J]sin2(^,), (4.9) 

^ ^ l<i<j<n i=l 

and the formula for /c-point correlation function is 

PuA^u ...M= det (i^2n+i(^., ^i))^j=i...... (4-10) 

The analogue of (4.4) reads 



K(^exp(t5^/(^,)) 

fe=i fe=i / 



(4.11) 



We will prove Theorem 2 for S0(2n). The proofs for S0(2n + 1) and Sp(n) 
are almost identical. 
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Proof of Theorem 2 The arguments from §1 imply that it is enough to 
prove 

Lemma 3 Let Q,n(/) be the £-th cumulant ofYTj=i fi^j)^'^ ^ 2. Then 



ki+...+ke=0 ^ 



+ G{-ki, . . . < const£ ^ 



(4.12) 

fci + . . . + fc^ = 

|fei|+...+|fc^|>n 

|A:i||/(A:i)| • . . . ■ |/>,)| + const^ J] |/>i)| • . . . ■ 

|fciH-...+|fe^|>n 

□ 

We start with the formula (2.6) which holds for general determinantal 
random point fields: 

(.i>i 

m 

{K2n-l{Xj, Xj^i) + K2n-l{Xj, —Xj^ij) clXi . . . dx^ 

(we always assume Xm+i = xi). 

E E- E 

m=l + ... + «„,= ^ ei=±le2=±l em=±l 

4>1 



(4.13) 



/ f^{x\) ■ ■■■■ f^^ixm) ■ Yl ^2n-i{xj, ej ■ Xj+i)dxi • . . . • dx^ 
Each term in the last sum with YliLi = 1 is equal to 

P mm 
/ f'{xi) • . . . • f^iXm) ■ n K2n-l{Xj,Xj+i) TT d{ei-i ■ Xi) = 

If ™ 

7^ • / f'{xi) ■ ■■■■ f"-{Xm) ■WK2n-l{Xj,Xj+i)dXi ■ . . . ■ dXr, 
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(we use the fact that f{x) is even). Combining these terms together we 
obtain the same expression as for | ■ Cg^2n-i{Ylf=L^ fi^j)) ^^e case of 
U{2n — 1), which gives vanishing contribution if ^ > 2. Finally we claim that 
the contribution from the terms with YYiL^ = — 1 can be bounded from 
above by 

const^ Yl \f{ki)\-...-\f{ke)\. 

\ki\+...+\ke\>n 

Indeed, the integral 

can be rewritten as 

^ n n 

Y ■■■ Y f^^{si- em- Sm) ■f^^{s2- ei- Si) ■ . . .■f^^{Sm- em-1- Sm-l) 



Consider the euclidian basis in and define fj = Cj — ej_iej_i, eo = 

Sm- The vectors form a basis in iff H^i — Then for any 

m-tuple {ti, . . . ,tm) there exists the only m-tuple (si, . . . , s^) such that 



tj = - ^j'l ■ J = 1, • • • We write f^^itj) = Yl f{hi+...+e,^,+i) ■ 

. . . ■ /(fc4+„,+£j, where the sum is over ki such that ki — tj. 

When we plug this into (4.13) we obtain a linear combination of 

/>i)....-/(U (4-14) 

It is easy to see that for |A;i| + . . . + < n the coefficient with the term 
(4.14) is equal to 



1 ^ 

L.y y 



1 P 

_ 



m ^l\■...■^m\ 

For + . . . + l/c^l > n the coefficient is bounded from above by some 
constant. This finished the proof of Lemma 3. □ 
Similar to §1 wc obtain the proof of Theorem 2 by applying the lemma 

to Y.U9iLn-{e,-e)). 
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